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SUMMARY 

General equations are developed for isentropic, frictionless, 
axi symmetric flow in rotating impellers with blade thickness taken 
into account and with blade forces eliminated in favor of the blade- 
surface function. It is shown that the total energy of the gas 
-relative to the rotating coordinate system is dependent on the 
stream function only, and that, if the flow upstream of the impeller 
is vortex- free, a velocity- potential function of only the radial 
and axial distances exists in the impeller. 

The characteristic equations for supersonic flow are developed 
and used to investigate flows in several configurations in order to 
ascertain the effect of variations of the boundary conditions on the 
internal flow and the work input. Conditions varied are prerotation 
of the gas, blade • burning rate, gas velocity at the blade tips, • 
blade thickness, and sweep of the leading edge. 


INTRODUCTION 

No method of practical importance has been described in the 
literature for computing the flow in the three-dimensional region 
between the blades of the supersonic compressor. This deficiency in 
design technique has been partly met by using two- dimensional theory 
to compute the variations of the relative flow from blade to blade, 
by assuming no radial flow components, and by neglecting spanwise 
variations. This procedure is probably satisfaqtory for designs with 
ratio of blade-root radius to blade-tip radius close to unity, and 
cylindrical hub and casing shape. When this ratio is small, the 
spanwise variations are significant. A method is described herein, 
which was devised at the NACA Lewis laboratory, for computing the 
magnitude of these variations, when blade-to-blade variations are 
neglected, by assuming a large number of blades in the impeller and 
dealing with the circumferential average flow. This technique, when 
combined with estimates of blade-to-blade variations, is useful for 
designs with closely spaced blades. 
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The usefulness of this approach may possibly be extended to super- 
sonic flow in a manner analogous to the method described in reference 1 
for subsonic flowj the axisymmetric solution is first found and then 
used as a base for computing the blade-to-blade variations. The par- 
ticular method used in the blade-to-blade calculation of reference 1 
is series expansion^ which is not applicable to supersonic flows 
because of the existence of solutions with discontinuous derivatives. 

In reference 2, the equations for internal flow with zero blade 
thickness are translated into the characteristic equations for the 
derivatives of the stream function, and it is shown there that the 
equations are hyperbolic if the relative velocity is supersonic, 
provided the ratio of the tangential to the axial velocity is a given 
field function. (This condition includes the special class of blades, 
containing radii from the axis of rotation) . If the tangential velocity 
is prescribed, then the meridional velocity component is the decisive 
factor in determining whether the equations are hyperbolic. 

The present report treats the velocity components directly after 
the manner of reference 3 for arbitrary preassigned blade shapes. The 
design problem is solved by finding the hub shape that will give a 
prescribed velocity distribution along the casing (or any other 
selected stream surface) . In the course of the solution, the tangen- 
tially averaged flow is completely determined in the interior of the 
impeller. 

The purpo'se of this paper, in addition to describing the method 
in a mathematical outline, is to examine the performance of the tech- 
nique in computed examples and to find if any peculiarities in the flow 
arise from the boundary Conditions imposed, and to develop a background 
of experience for the design of impellers of the type considered by exam- 
ination of the effect of systematic variations' of impeller design para- 
meters on the internal flow. The parameters varied are -inlet pre- 
rotation, blade turning rate, gas velocity distribution at the easing, 
blade thickness, and sweep of the blade leading edge. 

Approximately 160 hours of computing is required to work out about 
100 points for each of the impellers. This time can be reduced by 
application of the punched-card computing technique to part of the 
procedure. 

DEVELOPMEHT OF EQUATIONS FOR AXIALLY SYMMETRIC FLOW 

Equations describing the fluid flow are developed in convenient 
form by use of the equations of motion and continuity, the conditions 
of axial symmetry, and isentroplc state changes. The uniformly distri- 
buted blade forces, which are inconvenient in .calculations, are 
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eliminated from the equations by substituting the blade surface 
function; so that the condition of integrability is automatically 
satisfied. 


Equation of Continuity 
The continuity equation for steady flow is 

V '(pw) = 0 

(Symbols are defined in the appendix.) If axial symmetry is assumed, 
the equivalent in cylindrical coordinates is 

(pvr) + ^ (pwr) = 0 

TTip presence of blades serves to modify these equations if the space 
filled by the blades is taken into account. If the impeller contains 
N blades with thickness t measured in the direction of rotation, 
the ratio of volume available for fluid flow to the actual volume of 
any space element approaches (in the limit as the mmiber of blades 
increases indefinitely with Nt fixed) 

A ^ Et_ 
r ~ ” 2rcr 

The continuity equation is therefore obtained by reducing the specific- 
mass-flow terms pv and pw by the ratio A/r to accoiuit for the 
area reduction: 




(pwA) = 0 


( 1 ) 


Because none of the scalar quantities is a function of 0, equa- 
tion (l) implies the existence of a stream function \|r(r,z) such that 

= pwA 1 


= - pvA 


( 2 ) 


tg = 0 


(subscripts indicate partial derivatives) or 
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7\lf = pAj X V = pAj X W (2a) 


Equation of Motion 

Substitution of energy for pressijre. - The equation of motion for 
steady, nonviscous, isentropic flow in terms of the relative velocity is 

2 

V ^-W,XVXW + ivP = ttt^R + 2Wxn+F 

c p 

where co^R is the centrifugal force and 2W X Q the Coriolis force. 

For isentropic conditions, the substitution 

^VP=7h 

may be made, so that the equation is then 

„2 2 2 

w CO I* 

7 \ h + j = Wx(7XW + 2Q)+F = WX7XV + F (3) 

The quantity 


tt2 2 2 „2 

oor , V / V 

E = h + - — p — = h + - c 0 t(u + (ur) 


(4) 


is a measxire of the total energy of the gas particles in the relative 
system of coordinates and remains constant (as will be shown for 
frictionless flow) for any particle} energy is interchanged among the 
forms of enthalpy, kinetic energy, and potential energy ( - o?r^/2) 
resulting from the centrifugal- force field. In reference 2, E is 
shown to be constant for any particle for frictionless flow and no 
heat trsinsfer, with a finite number of blades in the impelfer. 

In order to show that E is constant for any particle in axially 
symmetric flow, the direction of the field force is first identified. 

If the blades are considered to increase indefinitely in nianber eind to 
retain their camber-line shape, the forces exerted by the blades are 
merely the pressure forces exerted normal to the blade surfaces, which 
in the limit approach uniformly distributed forces in the body of the 
fluid. The relative flow lines are always contained in one of the 
family of blade surfaces, so that the field or blade forces must be 
normal to the velocity. If the scalar product of equation (3) is taken 
with the relative velocity, the result is then 
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W • VE = 0 . 

which indicates that E does not change along any flow line. 

Elimination of blade force. - The blade force can be eliminated 
from the equation of motion if the equation for the blade surfaces is 
assumed known. The force may then be evaluated after the velocity com- 
ponents are known. 

Under the condition of axial symmetry, the equation for the family 
of blade surfaces 

a(r, 0, z) = constant 

may be simplified. If the constant is changed, there is a shift from 
one blade to another. Because of axial symmetry, this shift is equi- 
valent to rotating the original blade by a fixed angle. Therefore, the 
equation for the surfaces can be put into the form 

a = 0 + f(r,z) (5) 

Because the relative-velocity flow lines are .contained in the blade 
surfaces, W is normal to Va as well as normal to y^r. Therefore, 
the relation among the three vectors can be shown to be 

W = X Va (6) 

From the condition of axial symmetry, 

Se . _ 

30 = rj -VE = 0 

Because both VE and are normal to j and W, they must be 

parallel, and consequently E is known to be a function of f only 
(reference 4); that is, 

E = E(\lr) 

7E=wg 

Insertion of equations (4), (6), and (?) into equation (3) yields 
(without resolution of the term V x V) 
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= (^Vlr* V xvjva -(^Va. V X vjv^lr + F 

The blade pressure forces F are normal bo the blades and therefore 
parallel to vo-J "this relation thus includes only two nonparallel 
vectors Vo. and S 7 ^, so that the components must separately satisfy 
the relations 

F = - ^ V\jr • V X = - (r j • W X v ^ V)va- (8) 

and 


d\lr pA 


( 9 ) 


Special conditions. - Two special conditions of technical import- 
ance are considered: In hi^-speed machines^ blades with radial ele- 

ments are desirable in order to eliminate bending loads resulting from 
centrifugal forces. Another frequently occurring circumstance is that 
of zero vorticity upstream of the impeller. 


For blades with radial elements , the vector i lies in the blade 
surface and therefore 


i • va = f ^ = 0 

Therefore f is-a function of z alone. From equation (6), 


0 = W • Va = ^ + yfj. + wf2 


( 10 ) 


so that for radial blade elements 


u 


w r 


dz 


(10a) 


This relation provides considerable simplification of the equations 
developed later for solving the flow problem. 

For vortex-free flow upstream of the impeller, equation (3) is 
reduced for conditions there by setting n=F=^X V=0 and substi- 
tuting V for W to get 
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or 



= constant 


Also^ zero vorticity gives ^ for axial symmetry, 

„ _ S / X . / Sv Sw) k ^ / Z\ 

0 = 7 X V = - 1 ^ (u + mr) + j (ur + cor ) 

From the first and last terms and axial symmetry, 

p 

ur + cor = constant 

Then the value of the total energy E relative to the impeller is 


E = h + ^ 


(o^r^ 


y2 , 2\ 

= h + -g — co(ru + cor ) 


which is constant upstream of the impeller. Then hy equation (?) E 
remains constant in the impeller, so that the equation of motion (9) 
becomes 


V a • V X V = 0 


(11) 


If the system of coordinates used is r, z, and a for expressing 
equation (ll) , there results (reference 4) - 




(rV -7 a X k) 


r.a 


a 


(rV • i X 7a) 


z,a 


The subscripts indicate which variables are held constant. This equa- 
tion and the fact that the quantities in parentheses are functions of 


r and z only establish the existence of a potential function of 
r and z only ($(r, z)) such that 


a® 


z,a 


rV/7aX k 




lr,a 


rV • i X7 a 


= 0 

r,z 


Sa 
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The last relation indicates that with r and z constant ^ ^ is 
constant, and this is true regardless of the third independent 
variable, which may be taken as 6. The three equations for $ are 
then put into vector form: 

y ® = i(rV ’^a, X k) + k(rV * i X ya) ' (l2) 

Thp absolute velocity may then be shown to be 

V = y $ + r(u + cjur)ya (13) 

and from 

V • Va = W ’y a + air j • y a = CD 
equation (l3) gives for the moment of momentum p., 

„ =r(ut(m:) . ( 14 ) 

(yo-)^ 

SO that . 

■ (15a) 

The potential $ can be used as an alternative to the stream function 
in solving for the flow. Equation (13) indicates that the projection 
of V<I> on the blade surface is equal to the projection of the absolute 
velocity on the blade surface. The equation for \lr and that for $ 
are obtained in the following manner: The velocity V is eliminated 

from equation (9) by means of equation (6) and there is obtained 


/ r , \ o _n pA iE 

^ ^ j + ~ d^ 

which iB also applicable in the general case with nonzero vorticity 
upstream of the impeller. 
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In order to obtain an equation for *1*, equation (13a) is substituted 
into equation (2a) : 

Vt = pAiye X V = pAivQ xy$ + pAprye xyo- 

Then 


0 = yxv^ = VX [pAty 0 x ( y $ + ijya)] 

Although obtaining a solution in terms of \|r is feasible whether 
dE/dilf is zero or not, when dE/d\|r = 0 , ^ may prove more convenient. 
In regions where W^/a^ is nearly 1.0, because tw| is a double - 
valued function of the stream-function derivative p lw| , it will be 
uncertain whether the flow is subsonic or supersonic. A solution 
obtained from a network would req,uire a very fine spacing in regions 
where |W| = a. - 


General Equations in Cylindrical Coordinates 

In general, where f^ ^ 0 and dE/dilr ^ 0, the equations for the 
internal flow are solved simultaneously for the two velocity components 
V and, w, when u has been eliminated from the equations by using 
equation (10). The eqiiations are translated from their vector form 
to scalar notation by utilizing cylindrical coordinates. Then a 
further transformation is employed with characteristic quantities as 
the independent variables. 

The continuity equation (1) is first written in vector form (the 
condition of axial symmetry being used) and multiplied by r/pA: 

^V^p^w) = y.W + W.^ + |w. V(A/r) = 0 


For Isentropic flow 


VP _ Vh 
P a2 
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Introduction of E by equation (4) and use of W • V E = 0 produces 
„ W •yW^/2 . (s^rv . r( ^A/r hk/r\ ^ 

’ • * Ar ^ " -srj- ° 


a a 

This exi)ression finally reduces to 


- lz)''z - 2 ’-''z +"r) * - ^2)''r ‘ ^2 h\' h)^ * 




In cylindrical coordinates, equation (9) becomes 

^ ^ . dE 

^r - + ^r ^^z - ^^r = ^ 


(15) 


(16) 


If equation (l4) is now substituted into equations (l5) and (l6 ) y there 
results for continuity 


ruw f„ 

II - ^ t V, + 


rruv f^ /ruw f^ ^ 


^2 ru fj.v 
a2 ^ a2 


I (i + - i_ / + 

r a2 a2'^ 


^ P ^zz + ^rz + ^rr) 

I (A/r) + V ^ (A/r)] = 0 


and for the equation of motion 

(l + r^fz^] w^ + r^fzfr i^r ~ ^z) ~ (^' + ^^^‘r^) "^z + 


2„ 2 

r f 

z 


a /I eUw s / 


Sr \ f 


V ^z/j 


- 2 ^x^ - pA|S= 0 

r dijr 


(17) 


(18) 


cn '^'7 
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Use of Characteristic Variables 

Statement of basic equations. - The method of reference 3 is 
applied to equations (17) and (18) . The object is to find two alter- 
native equations with new variables | and t). Instead of z and r, 
such that each equation (corresponding to equations (l7) and (l8)) 
contains derivatives with respect to only one of the Independent 
variables. Two additional equations are required to find the position 
coordinates z and r in terms of the characteristic coordinates 
I and T|. There are two characteristic slopes dr/dz^ and 

the first corresponding to | varying and t) constant, the second 
to I constant and varying. If ^ represents either of these 
slopes and a represents either of the coordinates corresponding to the 
appropriate the results are 


= C Zr 


^ B 


■ Vz ■ ^ (\7x)2±/N/r2 (^)2 (m2 _ i)J 


(the upper sign of ± refers to the lower to and, 

Jw^j + GVjj + LZjj = 0 


where 


B = 1 + ^ ^2 


m 2 _ (u^ + + w^) 

a2 

(w - ruf^)^ 2 2x 

J H 2 (1 + r O 






(19) 

( 20 ) 

( 21 ) 


G = (m 2 - 1) ^ (^)2±^r2 (^)2 (m 2 . i) 

Sl 


>(21a) . 


D Sr2f 2 

d z 




dE 

+ cur)f^ - pA 
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Equation (20) for ^ shows that the system is hyperbolic if >1. 

In the homogeneous equations = D2 ** 0) , the coefficients 

and L_ (corresponding to and ^_) are zero; consequently, 

in the hodograph plane the slope of the characteristics would be 
given by ^ ~ “ 'g* In plane flow, J/G is a function of v and w 

only, and this equation may be integrated and the hodograph solu- 
tions tabulated in advance of application to any particular problem 
that may arise. In the present circumstances, however, J/G also 
involves the position coordinates, and ^ 0, Dg 0. A corre- 
spondingly easy method of solution by simple waves is therefore 
unavailable. Instead, the four (s = | or t]) , equations (19) and (21) 
are integrated simialtaneously in a step-by-step procediare. 

Method of computation. - In figure 1, conditions eire ‘ supposedly 


known at P^ 

and Pg, and conditions 

are to be computed at P^. 

Then 

equation (l9) 

in difference form is 




r3 - ri = (Z3 - 

2 l)^ + 

(22) 


"■3 - ''2 = (^3 ■ 


(22a) 


This system may be solved for r^ and z^. Then similarly for equa-' 
tion (21) , the difference form is 

J(w3 - w^) + G^ (V3 - v^) + (Z3 - Z]^) = 0 (23) 

J(w2 - Wg) + G_ (vg - Vg) + L_ (Z3 - Zg) = 0 (23a) 

The subscript + indicates that tj is constant and the subscript - 
indicates that I is constant. In equation (22) and (23a) the 
coefficients J, G, L are averages between values for the points 
and P3. For the first approximation, only the values at P^^ are 

known and these values are used. 

When approximate velocity components and position coordinates are 
thus determined for . P3, the coefficients J, G, L may be computed 

there and averaged for a better approximation of the position coordi- 
nates and the velocity components at P_. This process is repeated as 
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often as necessary to arrive at an answer of satisfactory accuracy. The 
smaller the interval, the quicker the convergence to the solution. In 
each stage of the process, if dE/dilr ^ 0, then must he computed 

from 

i|f3 - ilr = / + iif^)dz 

'^z 

From the value of' dE/dilr at may be found. 

• BOUIOARY COHDITIONS 
General 

The discussion in the previous section implies that boundary con- 
ditions are known on a curve C containing points such as and P^, 

and that the solution is then extended into a region containing points 
such as P^. The shape of the blades and E(\lr) are assumed known in , 

that region. On the curve C, there are known in addition the values 
for ^w and v, from which all other variables (such as dE/d\lf, \[f, u, 
p, a^, and so forth) may be computed. Suppose the initial curve is 
such that the characteristic curves intersect it at two points such as 
A and D or A and E in figure 2. Then by equation (21), there is a 
relation between the variables at A and D that does not permit the inde- 
pendent specification of the variables v and w at A and D. Simi- 
larly if two characteristics passing throiigh two boundary points reach 
a third boundary point, then another degree of freedom is lost. In 
figure 2, if two variables are prescribed in AB and EF, then one variable 
may be prescribed on BC and DE, and none on CD. 


Tip Conditions 

The supersonic impeller has a high pressure ratio per stage and a 
large air-handling capacity per unit of frontal area. In order to main- 
tain the large capacity, the tip radiiis at the entrance should be at 
least as large as the exit. In order to provide a maximum work output 
with a minimiam turning, the exit tip radius should be at least as large 
as the entrance tip radius. A reasonable shape for the tip streamline 
is therefore a surface of constant radius, although smaller radii between 
the entrance and the exit might 'be useful in order to affect the internal 
flow. If the casing shape is selected, then, as it is also a streamline. 
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the ratio v/w is specified. Selection of the velocity or a velocity 
component remains.. This selection may be made by consideration of 
the blade- to-blade velocity variations and by using some approximate 
method such as that outlined in the following paragraphs . 

In three-dimensional motion (finite number of blades), Kelvin's 
theorem states that if the absolute velocity is ifrotational upstream 
of the impeller, it will remain so inside the impeller; that is, 

VXV = 7XW + 2oik=0 

where the prime indicates velocities of the nonsymmetric flow. Then 


V t • 7 X V' = 0 


may be reduced to 

^ (¥' • q) = r + (u'+2cur) q • i (24) 

where Su'/Ss _is the arc-length partial derivative of u' along the 
direction of q. This equation is now integrated from the pressure 
s\irface of one blade to the suction surface of the next in the direction 
r = constant, z = constant. The variable 0 increases by the amount 
2jtA/Nr, so that 


(W • q)^ - (¥' . q)^ = 



+ 1 



u'de + 


2jtA 

IT 


If it is now assumed that the values of ^u'/ds and of u' averaged 
across the channel are equal ito the values obtained from the axisymmetric 
solution, and fiurther that -(¥' • q) = |¥' |, then 



if the relation 


(24a) 


V 1^1 = i • q 
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is utilized. In this equation all the velocity components without sub- 
scripts are obtained from the axisymmetric solution. If 

% 

|W’J + |Wp| = 2lwl 

is assumed or some other connection between , lWp|, jwj , and 

the impeller geometry, then the blade surface velocities at the casing 
may be obtained with that relation and equation (24a) . The distribu- 
tion of |w| at the casing may therefore be selected on the basis of 
obtaining a desired blade-surface velocity distribution or distribution 
of u. 

Suppose that, in some manner, all the variables have been selected 
at the case (curve BC in fig. 3(a)). The flow is then uniquely deter- 
mined in the region BP^C enclosed by the casing and two characteristic . 

lines. If the flow is also prescribed at the leading-edge surface BB*, 
then the flow is determined in the extended region B'BCP,P„. When the 

J_ ^ 

characteristic line BP^ falls inside this region (fig. 3(a)), no. diffi- 
culty results. If, however, the characteristic B’P^ intersects the 

casing inside the blade region (fig. 3(b)), then the velocity may not 
be prescribed on the cinrve BP^ (refer to discussion of fig. 2), or 

alternatively, only one of the velocity components may be prescribed 
on BB*. (BB' is not a streamline.) Similar considerations are valid on 
the curve CC at the exit. 

The selection of values for the variables v and w on the 
entrance edge BB' must be made by joining the solution for the flow in 
the region ABB 'A' with that for the flow in the blade region BCC'B' . 
Because the blades may deflect the air at BB' discontinuously and 
because the effect of thickness may be felt by' the air in a discontinu- 
oixs manner, the discontinuities in the flow should be evaluated. A 
method for making this estimate follows: 

The angle p is defined as the angle between the meridional 
projection of the blade leading edge and the z-axis (fig. 3(a)). The 
conditions outside the blade region are designated by subscript zero, 
inside the blade region by no subscript. Then across the surface BB' 
the equation of continuity yields 

Pq (^0 ^ ■ ^0 P) = fl - 1^1 p (w sin p - V cos p) 


(25) 
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The discontinuity in blade thickness can cause no force in the entrance 
surface BB', and the blades can exert no force in the blade surface. 
Therefore, no force and no discontinuity in momentum exist in the 
direction of the curves representing the intersection of these two 
siorfacesj that is 


(W - Wq) *(i cos p - k sin p) x 7a = 0 


or 


r(u - Uq) (fj. sin P + f^. cos p) = (v - v^) sin p + (w - w^) cos p 
There is also available 

u = - r ^wf 2 + vf 

If the process is assumed isentropic, then 

1 

1 


t ■ fer* ■ (‘ • 


2 2\^"^ 

V -‘V 
0 


2h, 


0 


is obtained by utilizing the fact that 


®“^0'^2" 2 “^■‘■2“ 2 


(26) 

(27) 


(28) 


and that 


There are then available equations (25), (26), (27), and (28) for 
finding p, u, v, and w. The initial values of fj. and f^ may be 

chosen from these equations to obtain desired relations between the 
velocity components outside and inside the blade regions. 


Applications of I^thod of Characteristics 

In the examples computed in this section, the conditions for which 
solutions are obtained vary. The purpose of these solutions is to find 
the effect of some variations of prescribed conditions on the resultant 
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flow. For all the examples, dE/d\lf = Oj the flow is vortex-free 
upstream of the impeller j at the leading and trailing edges the radial 
component v of the velocity is zero; and the blades contain radial 
elements (f^ = O) . Relative Mach number at the conjunction of the 

leading edge of the blade and the impeller casing is 1.94. Except in 
example VI, the value of the blade angle with respect to the axial 
direction is zero at the exit (f = 0) . The other conditions are 
summarized in the following table: 




Entrance 







Example 

1 

Pre- 

rotation, 

(u-wor)r 

Axial 

inlet 

velocity 

Tip 

speed 

Blade 

angle 

at 

tip 

(d.eg) 

Case 

velocity 

Turning 

rate 

Sweep- 

back 

(deg) 

Blade 

thickness 

I 

0 

0.85aj^ 

1.68a3_ 

-61.8 

(a) 

0 to 
maximxmi 
to 0 

0 

0 

II 

-0.3a^rt 

b.835aj^ 

1.44a^ 

-64.4 

1.94a^ 

0 to 
maximum 
to 0 

45 

0 

III 

-0.3a r 

X u 

0.835a^ 

1.44a^ 

-64.4 

1.94a^ 

Twice 

example 

II 

45 

0 

IV 


0.835a^ 

1.44a^ 

-64.4 

(^) 

Twice 

example 

II 

45 

0 

V 

-0.3a r, 
X "C 

0.835a^ 

1.44a^ 

-64.4 

(b) 

Twice 

example 

II 

45 . 

(c) 

VI 

-0.3a^rt 

0 . 835a^ 

1.44a^ 

-64.4 

(d) 

0 to 
maximTjm 

10 

0 


^etennined by conditions at leading edge to z/r^ = 0.356; held con- 
stant at 2.1a from z/r. = 0.356 to 1.0. 

•^Constant at 1.94a to = 0.25; decreased to 1.4a (= l.Ga^^) at exit. 

^Zero from z = 0 to = 0.25} from z/r^ = 0.25 to z/r^ = 0.5 

the correction for blade thickness is 


A = r exp 


' - — 1 - cos 8jt - 0.25] ' 

L " L l"t /JJ 


^Determined by leading edge; initial value, 1.94a. 
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Example I. - The results of the computation for the first example 
are summarized in figure 4, which shows streamlines, contours of con- 
stant Mach number and constant energy level, and some of the character- 
istics of the I -family, which intersect the leading edge. Because of 
the intersection of the characteristics from the leading edge with the 
tip of the blades, there is no control over the Mach number, as the 
blade shape was prescribed by 


d 0 df ■ 

- 3 — =-r-T— =-l, 
dz dz 


04 


(1 ^ cos , X j 


The Mach number at the case rises to a maximum of 2.1 at z/r^ = 0.356, 
and is then assigned a constant value. By equation (4), the stagnation 
enthalpy H (=h + V^/2), which measures the energy of the gas, 
increases with z at the rate 


Sh 


cur 


du 


At the case 


V = 0, and therefore, from equation (l4‘) and 
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Sh _ cur 
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This expression shows that the rate of work input is increased by 
increasing the rate of blade turning and decreased by increase in V. 
Fig\ire 4 indicates that the rate of tircning of the blade is too small 
in the beginning of the impeller and. that as a consequence nearly all 
the work input is concentrated in the back end of the impeller. This 
condition is aggravated by the acceleration of the gas in the initial 
portion (u|W|^<0 because u< 0) . Increasing the turning rate, 
however, will open the flow area and cause |W| to increase unless the 
area is controlled by blade thickness. 

A great disparity of work exists between the root and the tipj the 
values for AH/a^^ axe 0.65 and 2.5 at the exit. For isentropic com- 
pression, the corresponding ratios of outlet to inlet stagnation 
pressures are 2.48 and 9.0 for the root and the tip, respectively. The 
Mach number at the root at the exit (z = r.^., r = 0.507 r^) is in the 

neighborhood of 2.25, and a maximum value of about 2.5 is reached at 
z = r.^., r = 0.70 r.j.. Subsonic Mach numbers are reached in the region 

0<z/r^<0.3. The solution could not be extended into this region. 
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If the velocity were caused to decrease at the exit (z = r^) hy 
decreasing the velocity at the case (r = r^) in the rear part of the 
impeller (0.356 < z/r^ < 1.0) ^ then the mass-flow density at the exit 
would increase with the decrease in velocity and the required flow area 
would decrease. Termination of the huh would be ^xpected^ not at 
r = 0.507 r^, but at some larger, value where the work output would be 
higher. Another method of reducing the blade height at the exit and 
thereby reducing the work variation there would consist in relaxing 
the condition of v = 0 and utilizing some other relations that would 
correspond to lower Mach numbers. Two variables or relations could also 
be prescribed at the exit, but conditions on the blade tip (r = r,^.) 

could not be specified for some region near the back of the impeller. 

In this respect the entrance and exit conditions would be similar. 

Example II; Prerotation of gas. - In order to maintain supersonic 
flows throughout the impellers so that the flow coxald be computed, the 
relative entrance Mach number at the blade root for example II was 
increased by assuming an initial counterrotation. The rotative speed 
had to be reduced from 1.68a^ to 1.44a^ to maintain the work output, 

the relative Mach number at the leading-edge tip, approximately the same 
blade angle there, and the same axial velocity component at the inlet. 
This reduction had a further beneficial result of reducing the absolute 
Mach number at the exit for a given relative velocity axially directed. 
The condition of z = constant for the leading edge cannot be main- 
tained for shockless entry with radial blade elements and a counter- 
rotating free vortex entering the blade system. The sweep angle p 
chosen for this example was 45°, enough to permit free selection of 
W = 1.94a over the entire blade tip, because the characteristic curves 
from the leading edge do not interfere with boiindary conditions on the 
casing. The reduction in velocity at the casing reduces the blade 
height at the exit and the work-output variation there. 

In order to examine deviations from an example more nearly repre- 
sentative of a good design, the velocity at the blade tip for example II 
was not maintained the same as in example I, and consequently the differ- 
ences between examples I and II are not ascribable solely to inlet pre- 
rotation. Comparison of the effect on rotative speed (for the same 
turning at the tip and the same work output) and on diffusion at the 
blade ropts is, however, valid. 

The results of the computations are shown in figure 5. In this 
example, the turning rate is predetermined in the region -0.5< z/r,j. <0 

by the velocity distribution upstream of the leading edge, the condition 
of radial blade elements, and the slope of the leading edge. The 
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turning rate is continuous at z = Oj and satisfies the equation 

(0<z/r^ <l/3) 
(l/3<z/r^<l) 

Compared with the first example, the work-input rate at the blade tip 
is greater in the initial section; however, the work distribution is 
still unsatisfactory. The exit blade height is also much less, with 
the result that the work input is much more nearly equalized; at the 
root, the value of AH/a^ is 1.8, whereas, at the tip, the value is 
2.5. These values correspond to stagnation pressure ratios of 5.2 at 
the root and 8.35 at the tip. This variation could be reduced by 
sweeping back the trailing edge to provide space for additional work 
input at the hub streamline. This sweepback is more feasible with the 
smaller blade height and work variation of the second example than for 
the first example. 

Example III; Effect of turning rate. - In order to reduce the 
size of the impeller by increasing the turning rate, example III was 
computed with twice the blade turning rate of example II for z>0 
and with the same over- all work output. One result is the greater 
diffusion on the hub, where a minimum Mach number of 1.2 is reached 
instead of 1.4 as in the previous example (See fig. 6). The increased 
diffusion probably results from the increase in hub curvature, which 
in tiurn is caused by the shorter axial distances required for the 
changes. Another result is the increased height of the blade at the 
exit, and an increase in the exit Mach number from 2.3 to 2.7. The 
increased streamline curvature also accoiints for, the higher exit 
velocity and this higher velocity in turn explains the larger blade 
height. The characteristic line throxigh the casing at the exit 
(r = z = 0.5 r,j.) intersects the hub at z = 0.155 r^. The char- 
acteristics through the remaining portion of the hub (z >0.155 r^) 
are affected by conditions at the exit (z = 0.5 r,j.) where the veloc- 
ities increase with decreasing radius. The result is a large increase 

in velocity at the blade hub in the region z >0.155 r, . 

"C 

Example IV; Effect of diffusion at casing. - In fhe next example, 
all conditions remained the same as in the third example except that 
the prescribed velocity at the .case was varied smoothly from the 
value of 1.94a at z = 0.25 r,^. to 1.40a at z = 0.5 r,^.. The desired 
reduction has been obtained in blade height and work-output variation, 
but the discharge velocity has a large value of 2.6a (fig. 7). It may 


_d 

dz 


(l^) = 3.008 - 2.298 (3z - l)^ 
\dz / ' ' 

= 3.008 |l - 1/4 (3z - 1)^] 
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therefore be concluded that the alternative of directly prescribing a 
condition with low relative velocities at the exit would have been, 
more effective than prescribing conditions at the blade tip for 
reducing the blade height at the exit. 

A new phenomenon has also arisen in that the characteristics 
passing through z = 0.5 r^, r = 0.92 r^y and z = 0.5 r^, r = 0.893 r^ 

intersect at z = 0.397 r^., r = 0.832 r.^. At this point, the Maoh 

number changes discontinuously from 2.13 to 2.43 and the streamline 
inclination to the axis changes from 14.9° to 7.4°. This single point 
does not represent the full extent of this difficulty. Actually, there 
is a region bounded by a curve in which two solutions are obtained; 
each solution is determined by a different region of the boundary. Con- 
sequently, the boundary conditions are incompatable with obtaining a 
single-valued continuous solution. This situation is analogous to the 
Prandtl-Meyer flow, where the characteristics intersect at the corner. 

If a computation were started for a simple two-dimensional supersonic 
flow from a streamline shape and from the velocity distributed thereon, 
which were exactly the same as on one of the streamlines of a Prandtl- 
Meyer flow around a corner, then, when the solution was extended toward 
the corner, the characteristics would intersect at the corner in the 
manner of figure 7. The intersection on figure 7 can be taken to repre- 
sent a physical boundary to the fluid in that the streamline passing 
through that point wo\ild be taken for the hub shape. Only 90 percent 
of the flow (\lf = 0.9) could be accommodated in the impeller of 
example IV if continuous flow with the imposed boundary conditions were 
required. An obvious method of increasing the flow is to increase the 
slope of the characteristic lines .at the exit for r< 0.9 r . . This 

"G 

increase could be accomplished by decreasing the relative velocities 
there and discarding the condition v = 0 to obtain a shorter blade 
height . 

Example V; Increasing blade thickness. - In none of the examples 
computed so far was blade thickness considered. An obvious effect of 
increasing blade thickness is to decrease the velocity in a region 
where there is no independent control of the velocity (such as the 
region z<0.35 r^ in example I, fig. 4). This phenomenon provides 
the degree of freedom required in such regions to control the relative 
velocity |W| as well as the surface velocities |Wg| and |Wpl . In 

regions where the velocity |W| at the blade tips is under control 
(the entire tip sections for examples II, III, IV) , altogether dif- 
ferent effects may be exp)ected. The effect of thickness under these 
circumstances was investigated by computing another example with a 
cosine-like bump in the blade thickness starting at z = 0.25 r^. 
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(See the previously presented table of botindary conditions.) The maxi- 
mum thickness was reached at z = 0.375 r^. At r = r^ the thickness 

is 20 percent of the blade spacing^ and at the hub the value is 25 per- 
cent. The results of the computation are shown in figure 8. The effect 
of the bump is not felt upstream of the characteristic with positive 
slope passing through z = 0.25 r = r^ . The most obvious effect is 

a decrease of flow Mach number in the region of increasing thickness 
(and inward from the casing^ where the velocity is not preassigned), as 
might be expected from one-dimensional theory. At z = 0.25 r^ and 

z = 0.375 r_^, the change in blade height is about 10 percent as compared 

with example IV. Another effect of thickness is to be seen in the 
waviness of the hub streamline in the region 0.15< z/r^< 0.25 where 

the thickness has not yet started to increase. This bump may be traced 
in other streamlines along the characteristic rimning through the point 
of maximum thickness at the case (r = r^, z = 0.375 r^) • The effect 

of the bump on the hub is therefore to increase the local velocity by 
curvature in order to hold the required velocity at the. casing} this 
action from hub to casing takes place along the chairacteristic line as 
simple waves in two dimensions, and the velocity is seen to have 
increased all along this characteristic line. The intersection of the 
characteristics and the double-valued solution occur near the exit over 
a considerably larger portion of the flow than in the previous example. 
In this example also, if a lower velocity were prescribed at the exit, 
a unique solution would be obtained. It would also be feasible to carry 
out this idea and at the same time to prescribe another condition' (such 
as V = 0) provided that one were willing to lose control of the condi- 
tion at the back end of the casing (r = r^) in a manner analogous to 
that used at the front end of example I. 


Example VI: Effect of leading-edge sweep. - It is of more technical 

interest, however, to carry out this scheme at the entrance, for then 
the leading- edge sweep may be reduced from 50 percent of the axial depth 
of the machine to considerably less. A computation of the inlet flow, 
conditions for such an impeller is shown in figure 9, where the condi- 
tions at the leading edge are the same as for example III. The blade 


d^e 




shape is different, with the second derivative of angle 

dz"' dz"' 

constant. Also, the velocity at r = 1.0 cannot be prescribed. Only 
in the region enclosed by the casing and leading-edge characteristics 
was the computation completed. The result (fig. 9) shows a large 
increase in Mach number, especially at the blade tip. Because of 1}his 
increase, the blade turning accomplishes no work in the beginning of 
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the impeller, but shows a region of energy abstraction from the air; the 
initial negative rotational component becomes more negative as the air 
turns axially and the velocity increases as in example I. The situation 
is quite different at the back end wher? the blade angles approach zero 
and the tangential component does likewise; all the work input is con- 
centrated in the back end of this impeller as in the other examples . 
Postponement of work input to the back end of the impeller is to be 
expected in any case where the velocity is increasing as the blades 
turn toward the axial direction, or if the velocity decreases as the 
blades tiorn past the axial direction. Internal flow and work distribu- 
tion could be considerably improved in this example by increasing blade 
thickness in the region of large velocity increase. 


CONCLUSIONS 

General equations developed for isentropic, axisymmetric flows in 
rotating impellers with any given blade shapes showed that the blade- 
force term can be eliminated and solved for, if desired, after the flow 
solution is completed. Further, if the flow upstream of the impeller 
is vortex- free, then a potential function of only radial and axial 
distances exists, which may be used to find the projection of the abso- 
lute velocity on the blade surface. 

For supersonic flow, the method of characteristics was used to 
compute the flow for several examples with radial blade elements and 
vortex-free inflow. These examples show that 

(1) With zero radial inflow, the influence of the leading edge and 
the blade turning on the casing velocity in regions where this velocity 
cannot be prescribed in advance causes an increase in velocity aniS a 
reduction of blade loading and work input in the initial portion of the 
compressor . These effects ^an be compensated for to some extent by 
increasing blade thickness. 

(2) Sufficient leading-edge sweepback adds a degree of freedom to 
the boundary conditions at the blade tips, and therefore permits a more 
constant or a decreasing velocity inside the impeller. 

(3) Reduction in the exit velocity results in a shorter blade 
height and a more uniform distribution of work input at the exit. 

(4) Increase in the rate of turning causes more rapid hub curva- 
tures and consequently larger internal velocity gradients. Diffusion 
on the hub is therefore greater in the initial part of the impeller, 
and a lower velocity minimum is reached there, whereas the gas velocity 
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at the hub near the exit reaches larger vaj_ues. This situation requires 
larger blade heights at the exit with a larger exit velocity at the root 
and greater variation of work output from root to tip. Diffusion of the 
gas velocity at the tip will moderate this effect; as will also a relax- 
ation of the condition zero radial velocity at the exit and the substi- 
tution of a condition giving lower exit velocities. 

(5) Not every set of boundary conditions results in a continuous, 
single-valued solution in the region of interest. In the examples con- 
sidered, a relaxation of the condition zero radial velocity at the exit 
woTild have permitted another condition giving lower velocities there. 
This condition would have eliminated the region of- the double-valued 
solution and given a shorter blade and a more nearly equalized work 
output. 

(6) The effect of variations in blade thickness in regions where 
velocity is also prescribed at the blade tip is to change the velocity 
as in one- dimensional flow but not by the same amount, for there is 
some compensation by variations in blade height. The prescribed veloc- 
ity and thickness require hub curvatvires that cause velocity changes to 
be projected like simple waves on characteristic curves. It is there- 
fore possible to find regions of increased velocity when the blade 
thickness is increased. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 24, 1951. 
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APPENDIX - SYMBOLS 

The following symbols are used in this report; 

A r - Nt/2rt 

a sonic velocity of gas 

a^_ sonic velocity of gas at conjunction of leading edge and 

impeller case 

a* critical speed of gas upstream of impeller 

B ah abbreviation (see equations (21a)) 

Di^Dg abbreviations (see equations (21a)) 

E sum of stagnation enthalpy and potential energy of gas (specific) 

L w2 

in moving coordinates ^=h + - 

F distributed blade force per unit mass of fluid (a vector) 

f part of blade surface function dependent on r and z o nl y 

(a = 0 + f (r, z) ) 

G an abbreviation (see equations (21a)) 

H absolute stagnation enthalpy (h + V /2) 

AH work input to gas (rise in absolute stagnation enthalpy) 

h static enthalpy of gas (specific) 

i unit vector parallel to curve 0 =: constant; z = constant 

(i = 7r) 

J an abbreviation (see equations (21a)) 

j unit vector parallel to curve z = constant, r = constant 

(j = rV0) 

k unit vector parallel to curve r = constant, G = ^constant 

(k = 7z) . 
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an abbreviation (see equations (21a)) 
local relative ^ch number 

number of blades 
static pressure of gas 
unit vector parallel to W 
vector of r 
normal distance from axis of rotation 

distance from axis of rotation to blade tip at leading edge 

arc -length element parallel to W 

blade thickness measured in direction of rotation 

rotational component of relative gas velocity 

absolute gas velocity (a vector) [iv + j(u-KJor) + kv0 

radial component of gas velocity 

relative gas velocity (a vector) (iv + ju + kvr) 

axial component of gas velocity • 

distance measured parallel to axis of rotation 

function for blade surface (a(r,0;z) = constant) 

angle between meridional projection of leading edge and axis of 
rotation 

ratio of specific heats 

either or in general form of equations 

slope of curve for t v^iable, ri constant 
slope of cxarve for t] variable^ i constant 
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Ti one of characteristic coordinates for gas flow in impeller 

0 meridional angle (for cylindrical coordinates) 
p absolute moment of momentum (r(u + oar)) 

1 one of characteristic coordinates for gas flow in impeller 
p gas density 

a either rj or I in general form of equations 

$ velocity-potential function of r and z only that exists for 
initially irrotational absolute motion 

\lf stream function (y\|f = p AJxV) 

n vector form of oj 

0 ) angular velocity of impeller 

Subscripts : 

0 upstream of impeller 

p pressure surface of blade at tip 

s suction surface of blade at tip 

Subscripts r^ 0, z, a, q indicate partial differentiation. 
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Figure 2. - Relation between vShape of boundary and boundary conditions. 










Figure 5. - Flow characteristics of examp. 
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Figure 6. - Flow characteristics of example III. 
















NACA TN 2388 


35 



Figure 8. - Flow characteristics of example 
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